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In this work, employing the exchange-only orbital-dependent functional, we have obtained the 
optimized effective potential using the simple iterative method proposed by Kummel and Perdew 
[S. Kummel and J. P. Perdew, Phys. Rev. Lett. 90, 43004-f (2003)]. Using this method, we have 
solved the self-consistent Kohn-Sham equations for closed-shell simple metal clusters of Al, Li, Na, 
K, and Cs in the context of jellium model. The results are in good agreement with those obtained 
by the different method of Engel and Vosko [E. Engel and S. H. Vosko, Phys. Rev. B 50, 10498 
(1994)]. 
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I. INTRODUCTION 

In spite of the success of the local density approxi- 
mation (LDA)[ul and the generalized gradient approxi- 
mations (GGA)0, H for the exchange-correlation (XC) 
part of the total energy in the density functional theory 
(DFT) Q , it is observed that in some cases these approx- 
imations lead to qualitatively incorrect results. On the 
other hand, appropriate self-interaction corrected ver- 
sions of these ;ppr;ximations are observed^ to lead to 
correct behaviors. These observations motivate one to 
use functionals in which the self-interaction contribution 
is removed exactly. One of the functionals which sat- 
isfies this constraint is the exact exchange energy func- 
tional. Using the exact exchange functional leads to the 
correct asymptotic behavior of the Kohn-Sham (KS) po- 
tential as well as to correct results for the high-density 
limit in which the exchange energy is dominated. Given 
an orbital-dependent exchange functional, one should 
solve the optimized effective potential (OEP) integral 
equation!^ S 13 to obtain the local exchange potential 
which is used in the KS equations. Application of this 
integral equation to three dimensional systems 0, 0, ^ 
needs considerable technicalities and has some limita- 
tions. Recently, Kiimmel and Perdew 0, |^] proposed 
an iterative method which allows one to solve the OEP 
integral equation accurately and efficiently. 

In this work, using the exact-exchange OEP method, 
we have obtained the ground state properties of simple 
neutral A^-electron metal clusters of Al, Li, Na, K, and Cs 
with closed-shell configurations corresponding to N= 2, 
8, 18, 20, 34, and 40 (for Al, only iV = 18 corresponds to 
real Al cluster with 6 atoms). However, it is a well-known 



fact that the properties of alkali metals are dominantly 
determined by the delocalized valence electrons. In these 
metals, the Fermi wavelengths of the valence electrons 
are much larger than the metal lattice constants and the 
pseudopotentials of the ions do not significantly affect 
the electronic structure. This fact allows one to replace 
the discrete ionic structure by a homogeneous positive 
charge background which is called jellium model (JM). 
For closed-shell clusters^he spherical geometry is an ap- 
propriate assumption 0, 0, and therefore, we ap- 
ply the JM to metal clusters by replacing the ions of an 
A'^-atom cluster with a sphere of uniform positive charge 
density and radius R = {zNy^^Vs, where z is the valence 
of the atom and Ts is the bulk value of the Wigner-Seitz 
(WS) radius for valence electrons. For Al, Li, Na, K, and 
Cs we take rs=2.07, 3.28, 3.93, 4.96, and 5.63, respec- 
tively. 

The organization of this paper is as follows. In section 
nil we explain the calculational schemes. Section IIIII is 
devoted to the results of our calculations and finally, we 
conclude this work in section Hvl 



II. CALCULATIONAL SCHEMES 

In the JM, the total energy of a cluster with exact 
exchange is given by 

+ ^ / (t)i[n,n+];r)[n{r) - n+{r)], (1) 

in which 
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(r)</.*,(r')</>,.(r)</>,.(r') 



r — r' 
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and 



(l){[n,n+];r) = 2 dr 



,[n(r')-n+(r')] 



I r — r' I 

Here, the background charge density is given by 

3 

n-^-{r) = n9{R — r); n 



(3) 



(4) 



and n(r) is calculated from 



(5) 



where 0icr(r) are the KS orbitals obtained from the self- 
consistent solutions of the set of equations 



[h 



KSa 



.(r) = 0. 



In Eq.®, 



Veffair) = v{r) + VH{r) + Vxcrir) 



ri(r') 



(6) 



(7) 



(8) 



(9) 



All equations throughout this paper are expressed in Ry- 
dberg atomic units. 

To solve the KS equations, one should first calculate 
the local exchange potential from the exchange energy 
functional. This is done via the solution of the OEP 
integral equation. Recently, Kiimmel and Perdcw [r3.[T^ 
in a simple and elegant way have proved that the OEP 
integral equation is equivalent to the equation 



^^*,(r)(/),;,(r)+c.c. = 0, 

1=1 



(10) 



in which i/iio- are the self-consistent KS orbitals and -0^0- 
are orbital shifts which are obtained from the solution of 
the following inhomogeneous KS equations 



(11) 



with 



Sia- are the KS eigenvalues which satisfy Eq. , and in 
the right hand side of Eq. H12|l . VxaiJ^) are the optimized 
effective potential and 



2 



dr'^tkm^l, (13) 



= / dr(j)*^{r)vxa{r)(l)ia{r), 



(14) 



Uxza ^ J dr0*^(r)M„^(r)0i^(r). (15) 

At the starting point to solve the self-consistent OEP 
equations (|ll|l - (|15|l . the self-consistent KLI orbitals 
and eigenvalues are used as input. Then we solve Eq. 
(|ll|l to obtain the orbital shifts ipia- In the next step, we 
calculate the quantity 



(16) 



i=l 



the deviation of which from zero is a measure for the de- 
viation from the self-consistency of the OEP-KS orbitals. 
This quantity is used to construct a better exchange po- 
tential from 



{r) = vt'{r) + cSAr). 



(17) 



With this w^CT™ i^) and keeping the KS eigenvalues and 
orbitals fixed, we repeat the solution of the Eq. 1)11(1 . Re- 
peating the "cycle" ((TT|l . l(T7|) for several times, the 
maximum value of 5o-(r) will decrease to a desired small 
value (in our case down to 10^^ a. u.). After completing 
cycles, the f"^™ in conjunction with the KS orbitals are 
used to construct new effective potential to " iterate" the 
KS equations ©. The value of c in Eq. H17|l is taken 
to be 30 as suggested in Ref.0|. We have used 10 cy- 
cles between two successive iterations. These procedures 
are repeated until the self-consistent OEP potentials are 
obtained. 



III. RESULTS AND DISCUSSION 

Taking spherical geometry for the jellium background, 
and solution of self-consistent KS equations, we have ob- 
tained the ground state properties of closed-shell 2, 8, 
18, 20, 34, and 40-electron neutral clusters of Al, Li, Na, 
K, and Cs in the exact-exchange jellium model and com- 
pared the results with those of KLI and LSDA. 

To solve the KS and OEP equations for spherical ge- 
ometry we take 



and 



Y, 



(18) 



(19) 



Substitution of Eq. (O and Eq. (O into Eq. |(TT1) 
the inhomogeneous KS equation reduces to 



■^+e^a-veffa(r) 



kik + l) 



(20) 
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in which 



q^^{r)=qfXr) + €J{r) 



(2), 



(21) 



with 



and 



<lil^{'^) = [vxca{r) - VxciG + Uxcia] Xi<y{r), (22) 



,(2 



N„ h+lj 



3 = 1 l=\U-lj\ 



2/ + 1 



Xja{r)B„{i,j,l]r) 



x[I[ljmj,limi,lmj — nii)] . (23) 
The quantities B and / in Eq. 123(1 are defined as 



/I 



Ba{i,j,l]r)^ / dr'xia{r')Xjcr{r')—- 

Jr'=0 ^ + 

+ / dr'x^Ar')xMr')^ (24) 



(25) 

and the bar over imphes average over rrii and nij. 
Also, the expression for Uxia reduces to 



47r 



[I{ljmj, hrrii, Irrij — rui)] 



X / dr X',a{r)x3a{r)B^{iJ,l;r). (26) 
Jo 
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FIG. 1: Source terms gjlg o- ''■'^'^ 'i'i'=o o- atomic units for 
Na2. As is obvious, the two terms are equal and opposite in 
sign, so that the orbital shift for N = 2 vanishes and the KLI 
and OEP results coincide. 



In Fig^ the source term components q^^p ^ and q\^Q ^ 
are plotted as functions of radial coordinate. As is seen, 
they are equal and opposite in sign so that they lead to 
zero orbital shift, i.e., ^;=o,ct('') = 0. This result in turn 
leads to the coincidence of the KLI and OEP results. 

In Figs. Ela) and|2Ib) the self-consistent source terms 
qia{r) of Ea. H22|) are plotted as functions of radial coor- 
dinate for Nag and Liig, respectively. The corresponding 
orbital shifts are shown in Figs. OJa) and|3Ib). It 
should be noted that qi^ (r) and (r) must behave such 
that 



J dr Qia{r)(l)ta{r) 







and 



are satisfied. 
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FIG. 2: (a)-Source terms qia in atomic units for Nas, and (b)- 
for Lii8. For Nag, only / = and I — 1 orbitals are occupied 
for each spin component whereas, for Liig, the orbitals with 
Z = 0, 1, 2 are occupied. 

In order to solve the self-consistent OEP equations, we 
use the KLI self-consistent results as input. For the KLI 
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FIG. 3: Orbital shifts in atomic units (a)- for Nas and 
(b)- for Liis. 

calculations, we use [Eq.(23) of Ref.^3| with ipiai^) = 0]: 

vfaHr) = — — -^{0..(r)0*,(r)ti„,(r) 
2ncr[r) ^ 

^ ^ i— 1 

+ |0,;<j(r)|^(U:ricr - U,j,ia)) + C.C. (29) 

The self-consistent exchange potentials of Li2 and Alis 
are plotted in Figs. Efa) and^Jb), respectively. For com- 
parison, the LSDA exchange-correlation potentials are 
also included. One notes that in Li2 case, the KLI and 
OFF potentials are completely coincident whereas, in Al 
case, the KLI and OEP coincide only in the asymptotic 
region. On the other hand, the LSDA potential, because 
of wrong exponential asymptotic behavior, decays faster 
than the KLI or OFF, which have correct asymptotic be- 
haviors of 1/r. In the Al case, = 18 refers to the 
number of electrons which corresponds to the number 
n = 6 of Al atoms. 

In Figs. Ela) and |SJb), we have shown the self- 
consistent densities for Li2 and Alig, respectively. As 
in the potential case, for Li2 the KLI and OFF densities 
completely coincide whereas, in Alig the coincidence is 
only at the asymptotic region. 



FIG. 4: Exchange potentials in KLI and OEP and exchange- 
correlation potentials in LSDA, in Rydbergs, for (a)- Li2 and, 
(b)- for Ali8. Here 18 refers to the electrons which is equiva- 
lent to 6 Al atoms. In Li2 the KLI and OEP completely co- 
incide whereas, in Ali8 the coincidence occurs at asymptotic 
region. The LSDA has wrong exponential decay whereas, KLI 
and OEP have correct 1/r decays. 



In Table m we have listed the self-consistent calculated 
ground state properties of the closed-shell clusters of Al, 
Li, Na, K, Cs for N =2, 8, 18, 20, 34, and 40. For com- 
parison of our OFF results with those obtained by Engel 
and Vosko (EV)[l3, we have also included those results 
for Al, Na, and Cs. The EV results are based on gradi- 
ent expansion which, in principle, is valid only for slow 
variations of density as in a bulk solid. However, for fi- 
nite systems such as clusters or surfaces, the FV results 
may differ from the exact OFF results. Comparison of 
our OEP total energies with those of FV for Na clusters 
shows a difference of 0.002% on average. On the other 
hand, the FV exchange energies differ, on average, by 
0.001% and the average difference in eh is 0.08%. From 
the computational costs point of view, these quite small 
differences makes the FV method advantageous for cal- 
culations within above mentioned accuracies. 

Now we compare the total energies and the exchange 
energies in the KLI, OFF, and LSDA schemes. Compari- 
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FIG. 5: Densities in atomic units for Li2 and, (b)- for Alis. 
Here, as in potentials, we have full coincidence for Li2 and 
asymptotic coincidence for Alig. 



son of the total energies shows that the OEP energies, on 
the average, are 1.2% less than those of the KLI. We do 
not compare the total energies of OEP and LSD because 
in LSD there exist a correlation contribution. On the 
other hand, comparison of the exchange energies shows 
that on the average, the exchange energies in the OEP is 
0.33% more negative than that of the KLI whereas, it is 
9% more negative than the LSD. 

An other feature in OEP which should be noted is 
the contraction of the KS eigenvalue bands relative to 
those of KLI. The results in Tabled show that for all N, 



the relation A°^p < A^^^ holds. Here, A = eh - £l 
is the difference between the maximum occupied and 
minimum occupied KS eigenvalues. For N—2, we have 
A = 0. The results show that the maximum relative 
contraction, I A°^P - A^^^|/A^^^, is 2.6% which corre- 
sponds to Csis- 



IV. SUMMARY AND CONCLUSION 

In this work, we have considered the exchange-only 
jellium model in which we have used the exact orbital- 
dependent exchange functional. This model is applied for 
the closed-shell simple metal clusters of Al, Li, Na, K, and 
Cs. For the local exchange potential in the KS equation, 
we have solved the OEP integral equation by the itera- 
tive method proposed recently by Kiimmel and Perdew 
[T^ . By solving the self-consistent KS equations, we have 
obtained the ground state energies of the closed-shell 
clusters {N = 2,8,18,20,34,40) for the three schemes 
of LSD, KLI, and OEP. The KLI and OEP results are 
the same for neutral two-electron clusters. However, for 
N 2, the densities and potentials in the KLI and OEP 
coincide for large r values. The OEP exchange and effec- 
tive potentials shows correct behavior of l/r compared 
to the incorrect exponential behavior in the LSD. The 
total energies in the OEP are more negative than the 
KLI by 1.2% on the average. On the other hand, the 
exchange energies in the OEP is about 0.33% more neg- 
ative than that in the KLI whereas, it is about 9% more 
negative than that in the LSDA. The widths of the occu- 
pied bands, £h — £l in the OEP are contracted relative 
to those in the KLI by at most 2.6% which corresponds 
to Csig. In spite of the validity of the gradient expan- 
sion method for slow variations in density, comparison 
of our OEP results with those of EV shows an excellent 
agreement. 
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TABLE I: Absolute values of total and exchange energies as well as highest occupied and lowest occupied Kohn-Sham eigenvalues 
in Rydbergs. The LSDA total energies include the correlation energies as well. 
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1 


,8431 


3, 


,5490 


0, 


,3211 


0, 


2134 


34 


4, 


,7230 


5 


.9117 


0, 


,3030 


0, 


,1908 


3, 


1617 


6, 


,1552 


0, 


,3320 


0, 


2229 


3, 


,1662 


6, 


,1934 


0, 


,3295 


0, 


2214 


40 


5, 


,5338 


6 


.8879 


0, 


,2995 


0, 


,1701 


3, 


,6226 


7, 


,0565 


0, 


,3234 


0, 


1988 


3, 


,6247 


7, 


,0744 


0, 


,3230 


0, 


1994 



Cs 5.63 



2 


0, 


,2382 


0, 


,2875 


0, 


,1789 


0, 


1789 


0, 


,1907 


0, 


,3335 


0, 


,2669 


0, 


2669 


0, 


,1907 


0, 


,3335 


0, 


,2669 


0, 


,2669 


0, 


,1907 


0, 


,3335 


0, 


2669 


8 


1 


,0252 


1, 


1904 


0, 


,2271 


0, 


1816 


0, 


,7341 


1 


,2742 


0, 


,2833 


0. 


2376 


0, 


7345 


1 


,2778 


0, 


,2826 


0, 


,2378 


0, 


,7345 


1 


,2777 


0, 


2378 


18 


2, 


,3652 


2, 


7459 


0, 


,2490 


0, 


1768 


1 


,6290 


2, 


,8866 





,2846 


0. 


2144 


1 


,6307 


2, 


,9044 


0, 


,2823 


0, 


2139 


1 


,6307 


2, 


,9043 


0, 


2132 


20 


2, 


,6351 


3, 


,0268 


0, 


,2471 


0, 


1548 


1 


,7969 


3, 


,1446 


0, 


,2772 


0. 


1904 


1 


,7980 


3, 


,1553 


0, 


,2773 


0, 


1920 


1 


,7980 


3, 


,1553 


0, 


1925 


34 


4, 


,5646 


5, 


2538 


0, 


,2613 


0, 


1743 


3, 


,0932 


5, 


,4652 


0, 


,2851 


0. 


2007 


3, 


,0974 


5, 


,5020 


0, 


,2830 


0, 


1994 


3, 


,0974 


5 


,5020 


0, 


1974 


40 


5, 


,3452 


6, 


1206 


0, 


,2584 


0, 


1554 


3, 


,5445 


6, 


,2591 


0, 


,2770 


0. 


1787 


3, 


,5462 


6, 


,2788 


0, 


,2766 


0, 


1791 


3, 


,5465 


6, 


,2763 


0, 


1795 



"Data from Ref 18|. 

'Here, A'^=18 corresponds to Ale cluster and other N's do not 
correspond to a real Al clusters. 
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